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ABSTRACT 


An  experimental  apparatus  which  will  allow  the  study  of 
the  influence  of  ground  effects  on  a  heated  cylinder  placed 
in  a  crossflow  of  air  was  designed  and  constructed. 

The  effects  of  a  wall  on  the  heat  transfer  character- 
istics of  a  heated  cylinder  were  investigated  at  a  Reynolds 
number  of  153,000  for  cylinder  to  wall  spacings  ranging  from 
d/r  =  0.0  to  d/r  =  5.33  (maximum),  where  d  =  gap  width,  and 
r  =  cylinder  radius.   Pressure  coefficients  were  obtained  at 
five  degree  intervals  around  the  cylinder  for  each  experiment. 

Experimental  results  obtained  for  d/r  =  0.0  agreed  well 
with  previously  reported  investigations  of  a  cylinder 
attached  to  a  wall.   Results  obtained  for  d/r  =  5.33  agreed 
well  with  previous  work  done  with  free  cylinders.   Experiments 
conducted  at  intermediate  values  of  d/r  resulted  in  minimum 
and  maximum  average  Nusselt  numbers,  suggesting  that  gap 
spacing  may  exert  a  strong  influence  on  the  behavior  of  the 
flow  field. 
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I.   INTRODUCTIOI 


The  principal  objective  of  this  work  has  been  to  design, 
construct  and  test  an  experimental  apparatus  which  will  allow 
the  study  of  the  influence  of  ground  effects  on  a  heated 
cylinder  placed  in  a  crossflow  of  air.   Previous  work  of  this 
nature  has  concentrated  primarily  on  free  cylinders. 
Recently,  McComas  (1)  investigated  the  pressure  and  heat 
transfer  characteristics  of  a  heated  cylinder  attached  to  a 
wall.   McComas  determined  the  local  and  average  Nusselt 
number,  and  the  pressure  coefficient,  on  the  surface  of  the 
cylinder  for  Reynolds  numbers  varying  from  53,000  to  153,000. 
This  work  is  a  fellow -on  study  of  the  project  initiated  by 
McComas,  differing  in  that  the  apparatus  has  been  designed 
and  constructed  such  that  the  plane  surface  can  be  varied 
from  the  attached  position  to  approximately  3  cylinder 
diameters  away.   Due  to  the  limitations  of  the  wind  tunnel 
used  by  McComas,  (described  in  detail  in  (1)),  hotwire 
anemometry  and  flow  visualization  techniques  were  not  able 
to  be  used  in  his  investigations.   The  present  experimental 
apparatus  is  located  in  the  Mechanical  Engineering  Department 
wind  tunnel  in  Building  500  on  the  Naval  Postgraduate  School 
campus.   This  wind  tunnel,  described  in  detail  in  section  III, 
has  a  speed  range  of  10  to  175  miles  per  hour,  and  is  equipped 
with  a  traversing  mechanism  and  hotwire  anemometer  which  can 
be  used  to  obtain  velocity  information  throughout  the  test 


section.   The  tunnel  also  has  a  capability  for  flow  visual- 
ization using  smoke  filaments. 

Experiments  were  conducted  at  a  free  stream  Reynolds 
number  of  153,000  for  cylinder  toplate  spacings  of  d/r  =  0.0, 
0.25,  0.5,  1.0,  2.0,  4.0,  5.33  (maximum  separation)  where 
d  =  gap  width  and  r  =  culinder  radius  (Figure  1).   The  flow 
field  in  the  vicinity  of  the  cylinder  was  probed  using  the 
hot  wire  anemometer.   The  velocity  distribution  suggested 
that  flow  blockage  was  a  factor  in  the  experiment,  and  that 
the  Reynolds  number  should  be  based  on  a  corrected  free 
stream  velocity.   The  correction  produced  an  effective 
Reynolds  number  of  161,000.   Using  this  corrected  Reynolds 
number,  the  data  taken  at  d/r  =  0.0  agrees  with  the  results 
of  McComas,  and  the  data  at  d/r  =  5.33  agrees  with  the  work 
of  Field  (2). 

As  was  the  case  with  McComas,  experiments  with  the 
cylinder  attached  to  the  plate  showed  a  critical  flow  pattern 
to  exist  over  the  top  half  of  the  cylinder.   As  the  cylinder 
to  plate  spacing  was  increased,  this  flow  pattern  changed 
from  critical  to  subcritical  over  a  narrow  range  of  plate 
spacing  (at  about  d/r  =  0.09).   The  flow  over  the  bottom 
half  of  the  cylinder  exhibited  critical  flow  as  soon  as  a 
space  existed.   This  flow  also  changed  from  critical  to 
subcritical,  but  the  change  was  less  pronounced  and  it  was  at 
a  larger  cylinder  to  plate  spacing  (at  about  d/r  =  1.2). 
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Figure  1.   Schematic  diagram  showing  determination  of  d/r 
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Another  phenomenon  which  was  seen  to  occur  with  increasing 
d/r  was  the  gradual  shift  of  the  stagnation  and  laminar 
separation  points  over  a  range  of  15  to  20  degrees. 

In  calculating  and  plotting  the  average  Nusselt  number 
as  a  function  of  gap  spacing,  an  interesting  and  unexpected 
trend  was  revealed.   It  was  found  that  the  minimum  average 
Nusselt  number  did  not  occur  when  the  plate  was  attached 
to  the  cylinder,  but  at  a  d/r  =  0.25.   Further,  this  curve 
also  indicated  a  maximum  point  at  d/r  =  2.0.   Although  the 
current  data  are  incomplete,  (further  investigation  is 
needed  over  a  range  of  Reynolds  numbers)  the  Nu  versus  d/r 
curve  at  Re  =  153,000  suggests  that  gap  spacing  may  exert  a 
strong  influence  on  the  behavior  of  the  flow  field.   Future 
studies  employing  both  flow  visualization  as  well  as  heat 
transfer  measurements  would  aid  in  establishing  a  plausible 
explanation  for  the  trends  observed  in  the  curve  of  average 
Nusselt  numbers. 


12 


II .   BACKGROUND 

Although  a  great  deal  of  investigation  has  been  devoted 
to  the  flow  and  heat  transfer  characteristics  of  cylinders 
far  from  boundaries,  see  Meyer  (5)  and  Field  (2)  for  a  list 
of  references,  virtually  no  references  have  been  found 
relating  to  the  heat  transfer  and  flow  characteristics  of  a 
cylinder  placed  in  close  proximity  to  a  plane  surface.   The 
practical  applications  of  such  information  would  extend  from 
aiding  in  conducting  thermal  analysis  of  container  stored 
ordnance,  rockets  and  missiles  to  the  determination  of  the 
optimum  size  and  spacing  for  heat  exchanger  tubes. 

Although  McComas  (1)  experimentally  determined  local 
and  average  heat  transfer  coefficients  on  the  surface  of  a 
uniformly  heated  cylinder  attached  to  a  wall,  McComas  found 
no  references  to  either  theoretical  or  experimental  work  on 
this  subject.   McComas  proposed  a  theory,  based  on  theories 
that  had  been  postulated  for  forward  and  rear  facing  steps, 
that  quantitatively  explained  the  trends  in  his  heat  transfer 
data.   Other  than  McComas'  work,  no  further  references  have 
been  found  dealing  with  attached  cylinders. 

Since  the  fluid  flow  field  and  thermal  field  are 
intimately  coupled  in  a  forced  convection  process,  an  under- 
standing of  the  flow  field  is  mandatory  if  one  is  to  under- 
stand the  heat  transfer  process.   A  thorough  discussion  of 
the  fluid  flow  and  heat  transfer  characteristics  of  a  free 
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cylinder  can  be  found  in  (3),  and  a  proposed  but  as  yet 
unverified  description  of  the  flow  around  an  attached 
cylinder  can  be  found  in  (1). 
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III.   EXPERIMENTAL  APPARATUS 

A.    WIND  TUNNEL 

All  tests  were  conducted  in  a  subsonic,  open  circuit, 
wind  tunnel  manufactured  to  Navy  specifications  by  the 
Kenny  Engineering  Corp.,  of  Monrovia,  California  (Figure  2). 
The  overall  length  of  the  tunnel  was  40  feet,  with  an 
entrance  bell  5  feet,  5  inches  wide  by  7  feet,  5  inches  high. 
The  test  section  was  20  inches  wide  by  28  inches  high  at 
the  entrance,  expanding  to  31  inches  high  to  allow  for 
boundary  layer  growth.   These  dimensions  resulted  in  a 
contraction  ratio,  based  on  area,  of  10:1.   The  tunnel  was 
powered  by  a  75  HP,  1750  rpm  variable  speed  motor  driving  a 
45  inch  diameter  adjustable  pitch  vaneaxial  fan  capable  of 
supplying  air  at  10,000  to  70,000  cubic  feet  per  minute.   The 
tunnel  was  equipped  with  six  graded  mesh  screens  to  control 
turbulence  and  a  24  inch  manometer  tube  tapped  into  the 
stilling  chamber  and  test  section  for  velocity  measurement. 
The  tunnel  was  designed  to  operate  at  velocities  ranging 
from  10  mph  (14.67  ft/sec)  to  175  mph  (256.67  ft/sec). 
Turbulence  intensity  was  taken  using  a  hot  wire  anemometer 
and  found  to  be  approximately  0.2%  at  a  free  steam  velocity 
of  98  ft/sec  (see  Appendix  D). 
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B.    HEAT  TRANSFER  EXPERIMENT 

The  test  cylinder  was  essentially  three  cylinders 
attached  end  to  end  (see  Figure  3).   The  two  end  sections 
were  hollow  softwood  cylinders  supported  by  a  hardwood 
brace,  curved  to  match  the  cylinder  contour.   The  central 
cylinder  was  a  sheet  of  resistive  paper  (Temsheet)  fitted 
into  stepped  portions  of  the  end  cylinders,  and  attached  to 
the  end  cylinders  and  hardwood  brace  by  double  sided  adhesive 
tape.   The  three  cylinders  then  formed  one  cylinder  20  inches 
long,  to  correspond  to  the  test  section  width,  and  having  a 
diameter  of  3  inches,  chosen  to  allow  maximum  resolution  of 
the  liquid  crystals  without  excessive  blockage  within  the 
test  section.   The  assembled  test  cylinder  was  tightly 
packed  with  glass  wool  to  minimize  heat  transfer  into  the 
cylinder,  and  thereby  reduce  the  time  needed  to  reach  steady 
state  after  a  voltage  was  applied  to  the  paper. 

The  resistive  paper  used  was  a  commercially  available, 
carbon  impregnated  paper, Temsheet ,  with  a  nominal  thickness 
of  0.039  inches.   The  width  was  sized  to  conform  to  the 
stepped  circumference  of  the  end  cylinders  (9.422  inches)  and 
the  length  was  12  inches.   The  circumference  of  the  resistive 
paper  cylinder  was  marked  in  5  degree  increments,  in  two 
locations  2  inches  apart.   This  was  the  test  section  used  in 
all  heat  transfer  experiments.   The  entire  outer  surface  of 
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the  resistive  paper  was  covered  with  three  coats  of  NCR  S-43 
liquid  crystals  having  event  temperatures  of: 


red:  109. 8°F 
green:  111.0°F 
blue:    112. 3°F    (Field  (2)) 


Power  was  applied  to  the  paper  by  means  of  wires  lead 
inside  the  cylinder  through  one  end.   The  wires  were  attached 
to  aluminum  foil  tape  by  means  of  alligator  clips,  and  the 
edges  of  the  tape  were  painted  with  silver  paint  to  assure 
good  electrical  contact  with  the  paper.   The  nominal  distance 
between  the  tape/paint  was  4.8125  inches,  resulting  in  a 
surface  area  for  power  dissipation  of  45.34  in  . 

The  power  source  used  was  a  Lambda  regulated  power 
supply,  model  LK  345A  FM.   The  voltage  applied  by  the  Lambda 
was  measured  by  a  Fluke  8100A  digital  multimeter.   This 
multimeter  was  also  used  to  measure  the  resistance  of  the 
resistive  paper.   Air  stream  temperature  was  sensed  by  a 
copper-constantan  thermocouple  installed  in  the  stilling 
chamber.   The  thermocouple  output  was  read  on  a  Leeds  & 
Northrup  millivolt  potentiometer. 

The  "wall"  used  in  these  experiments  was  manufactured 
from  1/2  inch  plexiglass.   This  platewas46  inches  long,  the 
maximum  practical  length,  to  minimize  the  influence  of  the 
trailing  edge  on  the  flow  patterns  around  the  cylinder.   The 
width  of  the  plate  was  20  inches,  to  conform  to  the  test 
section  width,  and  it  was  beveled  on  the  under  side  of  the 
leading  edge  at  an  angle  of  approximately  12  degrees. 
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The  plate  and  cylinder  were  placed  in  the  tunnel 
horizontally,  and  in  such  a  way  that  the  leading  edge  of  the 
plate  extended  four  cylinder  diameters  (12  inches)  ahead  of 
the  cylinder  centerline.   The  cylinder  was  attached  to  the 
test  section  walls  by  portions  of  the  end  cylinders  pro- 
truding through  holes  cut  for  that  purpose.   The  plate  was 
capable  of  being  moved  vertically  by  means  of  a  threaded 
rod  extending  through  the  tunnel  floor,  and  a  rear  brace 
which  locked  in  position  by  means  of  a  split  nut  (Figure  4). 
The  plate  was  capable  of  being  located  at  any  position 
relative  to  the  cylinder,  from  attached,  to  a  maximum 
possible  distance  of  8  inches  away  (Figure  1). 

C.    PRESSURE  COEFFICIENT  EXPERIMENT 

The  experiments  conducted  to  find  the  variation  of 
pressure  around  the  circumference  of  the  cylinder  used  the 
previously  described  plexiglass  plate  as  the  wall  and  a 
3  inch  diameter  aluminum  cylinder  fitted  with  four  pressure 
taps  90  degrees  apart  (Figure  5).   The  circumference  of  the 
cylinder  was  scribed  in  5  degree  increments.   Plastic  tubing 
led  from  the  taps  through  one  end  of  the  cylinder  to  a 
manometer  bank. 

The  manometer  bank  consisted  of  several  tubes  that  were 
inclined  at  an  angle  of  30  degrees  from  the  vertical.   All 
four  cylinder  pressure  taps  as  well  as  the  stilling  chamber 
and  test  section  pressures  were  read  from  this  manometer 
bank . 
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IV.   EXPERIMENTAL  PROCEDURE 

The  present  work  was  conducted  utilizing  different  test 

facilities  than  those  used  in  previous  investigations.   The 

wind  tunnel  was  smaller,  and  of  the  open  circuit  variety. 

Due  to  experimental  considerations,  the  cylinder  size  and 

orientation  were  different,  and  the  problem  of  varying  the 

relative  separation  of  the  plate  and  cylinder  was  unique  to 

this  work.   In  order  to  assure  that  meaningful  and  accurate 

data  would  be  obtained,  it  was  necessary  to  compare  initial 

results  with  previous  work.   Since  both  Field  and  McComas 

performed  experiments  at  a  Reynolds  Number  of  153,000,  it  was 

decided  to  check  the  present  apparatus  at  that  Reynolds 

Number.   For  this  reason,  an  air  stream  speed  of  98  ft/sec 

was  used  for  comparison  experiments.   Because  of  time 

limitations,  all  of  the  data  for  the  present  work  were  taken 

at  98  ft/sec,  corresponding  to  a  Reynolds  Number  of  153,000, 

where  Reynolds  number  is  defined  as: 

U   d 
Re   =   

Re 

U 


Reynolds  number 

free  stream  velocity  at  test  section 
entrance,  ft/sec 


cylinder  diameter,  ft 

2 

kinematic  viscosity  of  air,  ft  /sec 
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A.    HEAT  TRANSFER  EXPERIMENT 

The  same  experimental  procedure  used  by  Field  and 
McComas  was  used  in  the  present  experiments.   For  complete- 
ness, a  description  of  the  procedures  will  be  repeated  here. 

Prior  to  the  start  of  the  experiment  a  sufficient 
voltage  was  applied  to  the  cylinder  to  cause  the  liquid 
crystals  to  reach  their  event  temperature  and  change  color. 
This  voltage  was  usually  from  12  to  15  volts.   After  about 
a  half  hour,  time  to  allow  the  cylinder  and  glass  wool 
packing  to  come  to  a  steady  state  temperature,  the  power 
supply  was  turned  off  and  the  multimeter  was  switched  from 
DC  voltmeter  to  ohmmeter.  In  this  manner  an  accurate  value 
of  the  paper  resistance  was  determined  at  a  temperature  which 
corresponded  to  the  liquid  crystal  event  temperature.   The 
multimeter  was  then  switched  back  to  DC  voltmeter  and  the 
wind  tunnel  was  started  and  brought  up  to  speed. 

The  wind  tunnel  velocity  range  and  test  section  velocity 
profile  were  examined  prior  to  obtaining  experimental  heat 
transfer  data  (see  Appendix  D).   It  was  determined  that  the 
velocity  profile  was  uniform  at  the  test  section  entrance, 
and  that  the  turbulence  intensity  was  approximately  0.2%. 
After  numerous  checks  with  the  hot  wire  anemometer,  it  was 
decided  to  use  the  installed  wall  manometer  to  set  the  tunnel 
speed  of  98  ft/sec. 
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After  the  tunnel  reached  operating  speed,  the  voltage 
was  adjusted  such  that  the  liquid  crystal  color  bands  would 
appear  at  various  angular  positions  around  the  cylinder. 
The  first  color  band  to  appear  would  be  at  the  lowest  local 
Nusselt  number  on  the  cylinder.   For  tests  with  a  d/r  greater 
then  0.25  (where  d  =  distance  from  plate  to  cylinder  rim,  and 
r  =  cylinder  radius)  this  was  always  the  laminar  separation 
point  (Figures  13-19).   As  voltage  was  increased,  the  color 
lines  would  appear  in  one,  or  as  many  as  5  or  6,  angular 
locations  all  having  equal  Nusselt  numbers.   In  this  manner 
the  cylinder  surface  temperature  could  be  compared  to  the 
air  stream  temperature  as  determined  by  a  thermocouple 
mounted  in  the  wind  tunnel  stilling  chamber.   Using  this 
temperature  difference,  and  calculating  the  power  dissipated 
from  the  voltage  applied,  the  local  heat  transfer  coefficient 
could  be  calculated  for  the  circumference  of  the  cylinder. 
Because  of  the  way  the  cylinder  was  constructed,  readings 
within  15  degrees  of  the  seam  were  considered  unreliable. 
The  majority  of  Nusselt  number  data  were  collected  with  the 
seam  located  at  180  degrees.   The  Nusselt  number  in  the 
vicinity  of  180  degrees  was  determined  by  turning  the  cylinder 
such  that  the  seam  was  facing  forward  at  0  degrees.   Several 
experiments  were  also  made  with  the  seam  at  positions  of 
135,  225,  and  315  degrees.   The  final  position  of  0  degrees 
was  chosen  for  the  seam  after  it  was  observed  that  the  value 
and  angular  location  of  the  Nusselt  number  at  the  laminar 
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separation  point,  and  in  the  vicinity  of  130  and  210  degrees, 
was  consistent  with  the  values  obtained  with  the  seam  at 
180  degrees. 

For  all  experiments,  distinct,  repeatable  values  of  the 
Nusselt  number  versus  angular  location  were  obtained.   The 
gradient  as  determined  by  the  color  bands,  was  sharp,  with 
the  exception  of  a  region  near  the  forward  stagnation  point, 
which  appeared  nearly  isothermal.   The  gradient  was  also 
indistinct  in  the  turbulent  boundary  layer,  where  the 
"scrubbing"  action  of  cooler  air  caused  the  Nusselt  number 
to  change  with  time.   Therefore,  the  Nusselt  number  shown 
for  the  turbulent  regions  must  be  considered  to  be  an  average 
val ue. 

B.    PRESSURE  COEFFICIENT  EXPERIMENT 

With  the  pressure  cylinder  mounted  in  the  wind  tunnel, 
the  manometer  bank  was  leveled,  pressure  leads  checked  for 
tightness,  and  the  wind  tunnel  brought  up  to  a  speed  of 
98  ft/sec,  corresponding  to  a  Reynolds  number  of  153,000. 
The  cylinder  was  positioned  such  that  the  pressure  taps  were 
initially  at  0,  90,  180,  and  270  degrees.   The  manometer 
readings  for  the  four  taps,  the  stilling  chamber,  and  the 
test  section  were  recorded.   The  cylinder  was  rotated  in 
5  degree  increments  until  pressure  readings  were  obtained  for 
the  entire  circumference. 
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The  pressure  coefficient  was  determined  as  follows 
P-P      P-P 


hpU' 


P  -P 
s 


where  : 


pressure  coefficient 

density  of  air 

air  stream  velocity  (test  section  entrance) 

pressure  on  cylinder  surface 

stagnation  pressure  (stilling  chamber) 

air  stream  static  pressure  (test  section 
entrance ) 
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V.   RESULTS 

Initial  heat  transfer  and  pressure  coefficient  experi- 
ments were  made  with  the  plate  and  cylinder  attached 
(d/r  =  0.0).   A  thin  felt  strip  was  placed  between  the  two 
to  assure  a  tight  seal.   The  data  obtained  in  this  configura- 
tion agreed  well  with  the  results  of  McComas  (see  Figure  6). 
It  can  be  seen  that  the  trends  in  the  data  are  identical, 
that  local  maximum  and  minimum  values  are  within  experimental 
uncertainty  (see  Appendix  B  for  uncertainty  analysis),  and 
that  a  local  maximum  and  local  minimum  were  again  observed 
near  the  point  of  attachment  of  the  plate  and  cylinder.   This 
is  further  evidence  of  the  trapped  vortex  phenomenon  as 
proposed  by  McComas. 

The  next  series  of  tests  were  made  with  the  plate  as 
far  removed  from  the  cylinder  as  possible,  d/r  =  5.33,  or 
d  =  8  inches.   As  there  appeared  to  be  some  blockage,  the 
method  of  Pope  (4),  as  used  by  Meyer  (5),  was  used  to 
determine  the  necessary  correction: 


Re  _   =(1.0+  e)Re 
corr    v       ' 


where : 


e  =  e  .  +  e  . 
sb    wb 


and  e  ,  is  the  solid  blockage  correction  factor 

esb    3  V 
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and  e  .  is  the  wake  blockage  correction  factor 

WD 


e  .  =  1/4  (D/H)  Cn 

WD      '    v  '  '        D 


where: 


D  =  cylinder  diameter,  inches 

r  =  cylinder  radius,  inches 

H  =  total  tunnel  height,  inches 

Cn  =  drag  coefficient  (see  Meyer  p.  107) 


Therefore : 


e  ■  |  0-£)2   +  \   (fy)  (1.2) 


=  0.05 


and 


Recorr  =  (1'05)  153>000 


Re^v>     =      161   >000 

corr 


Using  the  Re„_„  to  compare  the  results  of  this  experi- 
a       corr       r  r 

ment  with  the  work  of  Field,  it  was  found  that  the  results 
compare  favorably  with  Field  and  with  the  theory  of  Schuh  (6) 
(see  Figure  7).   It  should  be  noted  that  Field  experienced 
critical  flow  at  Re  =  149,000,  while  the  experiments  made  in 
this  work  were  exclusively  subcritical  for  d/r  =  5.33.   This 
phenomenon  could  be  attributed  to  higher  free  stream  turbulence 
or  cylinder  deformation,  either  of  which  could  have  induced 
critical  flow  in  Field's  experiments. 
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Figure  7.   Comparison  of  the  heat  transfer  results 

for  d/r  =  5.33  and  Re     =  161,000  with 

corr 

the  results  of  Field  at  Re  =  148,000 
and  with  the  theory  of  Schuh 
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Figures  13  through  19  are  polar  plots  of  the  local 
Nusselt  number  and  pressure  coefficient  for  values  of  d/r 
ranging  from  0.0  to  5.33.   In  examining  the  seven  plots  for 
significant  trends  or  unusual  phenomena,  two  items  of 
interest  were  immediately  apparent:  one  was  the  change  from 
critical  to  subcritical  flow,  the  other  was  a  shift  in  the 
location  of  the  stagnation  point  and  the  laminar  separation 
points  as  the  plate  was  moved  away  from  the  cylinder 
(see  Figure  8) . 

The  transition  from  critical  to  subcritical  flow  was 
investigated  in  order  to  ascertain  the  d/r  ratio  associated 
with  the  transition.   As  the  flow  was  critical  with  the 
plate  attached,  and  subcritical  over  the  top  at  d/r  =  0.25, 
the  first  transition  point  of  interest  was  located  in  the 
interval  0.0  <_   d/r  <_  0.25.   By  alternately  moving  the  plate 
and  adjusting  the  voltage,  the  transition  point  from  critical 
to  subcritical  flow  over  the  top  of  the  cylinder  was  found 
to  be  approximately  0.08  <_  d/r  <_   0.10. 

As  the  flow  over  the  bottom  of  the  cylinder  underwent 
a  transition  from  critical  to  subcritical  flow  in  the  interval 
1.0  <  d/r  <  1.5,  this  region  was  investigated  further.   Using 
the  same  procedure  of  carefully  adjusting  spacing  and  voltage, 
the  transition  was  found  to  occur  at  approximately 
1.17  <  d/r  <  1 .25. 
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It  should  be  noted  that  while  both  flows  described 
above  are  definitely  critical,  there  is  significant 
difference  in  the  degree  of  criticality  on  the  top  and 
bottom  of  the  cylinder.   This  can  be  seen  from  the  magnitudes 
of  the  Nu  on  the  data  shown  in  Figures  9  and  10. 

In  order  to  glean  further  information  from  the  data, 
the  average  Nusselt  number,  Nu,  was  computed  from  the 
equation: 


2tt 


Nu  = 


-   I  Nu 

TT    J 


de 


which  was  numerically  integrated  by  computer  using  Simpson's 
Rule  of  Quadrature  with  a  step  size  of  5  degrees.   The 
results  are  shown  in  Figure  11a. 

It  is  interesting  to  note  the  unexpected  shape  of  the 
curve.   Nu~  drops  to  a  minimum  value,  then  increases  to  a 
maximum,  and  then  gradually  decreases  to  a  nearly  constant 
value  as  d/r  approaches  relative  infinity.   These  results 
prompted  further  investigation. 

In  order  to  help  determine  the  influence  of  the  local 

Nu  on  Nu"  a  plot  was  made  of  Nu     for  the  forward  (-90  to 
_     r  max 

90  degrees)  and  rear  (90  to  -90  degrees)  portions  of  the 
cylinder.   The  results  are  shown  in  Figure  lib  and  exhibit 
the  same  general  shape  as  Nu".   The  drop  in  Nu  from  d/r  =  0.0 
to  0.25  is  explained  from  this  plot  when  the  sharp  drop  in 
Nu     on  the  front  of  the  cylinder  from  d/r  =  0.0  to  0.25  is 
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d/r  -  0.6    x 

d/r  «  0.25   . 


Figure  9.   Comparison  of  the  local  Nusselt  number  at 
d/r  0.0  with  the  local  Nusselt  number  at 
d/r  =  0.25. 
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d/r  =0.0     x 
d/r  =0.25   * 


Figure  10.   Comparison  of  the  pressure  coefficient  at 
d/r  -   0.0  with  the  pressure  coefficient  at 
d/r  =  0.25. 
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Figure  11a.   Average  Nusselt  number  versus  d/r 
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Figure  lib.   Maximum  Nusselt  numbers  on  the  front  and 
and  back  of  the  cylinder  versus  d/r. 
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noted.   (The  maximum  local  Nil  for  d/r  =  0.0  occurs  in  the 

separation  bubble,  in  the  vicinity  of  90  degrees,  and  is  not 

used  in  the  plot  of  Nu    .  ) 

Further  investigation  of  the  Nu  curve  lead  to  an  attempt 

to  see  how  the  voltage  applied  to  the  cylinder  would  vary  if 

one  location  on  the  cylinder  was  maintained  at  a  constant 

temperature  as  the  plate  spacing  was  varied.   The  locations 

150  to  165  degrees  were  chosen  for  this  experiment.   For  each 

d/r  the  voltage  was  adjusted  until  a  distinct  color  band 

appeared  at  these  locations.   Plotting  the  voltage  applied 

versus  d/r  resulted  in  a  curve  of  similar  shape  to  Nu  and 

Num^„  (see  Figure  12). 
max        3 

In  order  to  verify  the  curvature  of  the  Nu"  plot,  a 
complete  set  of  data  was  taken  at  d/r  =  3.0.   This  data  is 
not  presented  on  a  separate  polar  plot  but  NTT  was  calculated 
and  the  value  obtained  enabled  the  shape  of  the  Nu  curve  to 
be  verified. 

A  possible  explanation  for  the  observed  variation  in 
Nu"  could  be  related  to  the  amount  of  "cold"  free  stream  air 
which  has  access  to  the  cylinder  surface.   At  d/r  =  0.0,  the 
flow  is  critical  due  to  the  fact  that  all  of  the  air  approach- 
ing the  cylinder  is  forced  to  go  "up  and  over."   As  soon  as 
the  plate  is  separated  the  flow  has  two  paths.   As  soon  as 
space  is  sufficient  to  provide  a  favorable  path  without  large 
pressure  gradients,  enough  of  the  air  approaching  the  cylinder 
passes  beneath  to  cause  the  flow  over  the  top  to  transition 
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to  subcri tical .   Now  the  air  passing  below  the  cylinder  is 
experiencing  critical  flow,  possibly  due  to  the  nozzle 
effect  produced  by  the  rapidly  decreasing  flow  area  below  the 
cylinder.   The  heat  transfer  in  this  region  does  not  increase 
due  to  the  fact  that  much  air  has  already  been  heated  passing 
over  the  front  of  the  cylinder  (i.e.,  the  thermal  and 
hydrodynamic  boundary  layers  are  a  significant  fraction  of 
the  gap  size). 

As  d/r  is  increased  beyond  0.25,  the  larger  opening 
allows  more  cold  air,  outside  the  boundary  layer,  to  approach 
the  bottom  of  the  cylinder.   This  cold  air  carries  off  some 
heat  by  conduction  and  also  mixes  with  the  warm  air  in  the 
separation  bubble  and  turbulent  boundary  layer,  resulting  in 
increasi  ng  TTu . 

The  gradual  decrease  in  Nu  from,  a  high  value  at 

d/r  =  2.0  is  not  understood  at  this  time.   It  was  expected 

that  Nu  would  approach  some  constant  value  as  the  influence 

of  the  plate  was  reduced.   It  can  be  seen  from  Figure  11a 

that  the  uncertainty  band  nearly  encompasses  the  decrease 

in  Nu.   Although  the  behavior  of  the  Nu    and  voltage  curves 
3  max  3 

suggest  that  the  trend  may  be  genuine,  further  investigation 

is  needed  to  determine  if  Nu  actually  has  a  peak,  or  in  fact 

approaches  a  constant  value  from  some  minimum  point  at  a 
small  d/r. 
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Figure  13.   Polar  plot  of  Nusselt  number  and  pressure 
coefficient  versus  angle  for  d/r  =  0.0. 
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Figure  14.   Polar  plot  of  Nusselt  number  and  pressure 
coefficient  versus  angle  for  d/r  -  0.25. 
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Figure  15.   Polar  plot  of  Nusselt  number  and  pressure 
coefficient  versus  angle  for  d/r  =  0.5. 
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Figure  16. 


Polar  plot  of  Nusselt  number  and  pressure 
coefficient  versus  angle  for  d/r  =  1.0. 
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Figure  17. 


Polar  plot  of  Nusselt  number  and  pressure 
coefficient  versus  angle  for  d/r  =  2.0. 
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Figure  18.   Polar  plot  of  Nusselt  number  and  pressure 
coefficient  versus  angle  for  d/r  =  4.0. 
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Figure    19. 


Polar  plot  of  Nusselt  number  and  pressure 
coefficient  versus  angle  for  d/r  =  5.33. 
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V I .   CONCLUSIONS  AND  RECOMMENDATIONS 

The  primary  objective  of  this  work,  to  design,  construct 

and  test  an  experimental  apparatus  which  will  allow  the  study 

of  ground  effects  on  a  heated  cylinder  placed  in  a  cross  flow 

of  air,  has  been  achieved.   The  comparisons  made  with 

previous  work  are  within  experimental  uncertainty,  and 

indicate  that  the  results  of  further  experiments  made  with 

this  apparatus  can  be  relied  upon.   The  liquid  crystal 

technique  for  heat  transfer  analysis  has  proved  to  be  an 

excellent  method  for  determining  quantitative  heat  transfer 

data.   The  results  of  the  one  plate  geometry  and  Reynolds 

number  investigated  in  this  work  are  shown  in  the  plots  of 

Nu  and  C   versus  theta.   Figure  8  shows  the  smooth  shift 
P 

in  the  location  of  key  points  which  characterize  flow 
phenomena  (stagnation  point,  and  laminar  separation  points) 
and  makes  it  possible  to  approximate  the  value  of  Nu  for 
configurations  other  than  those  investigated. 

Future  work  in  this  area  should  include  a  further 
investigation  of  the  variation  of  NU  with  increasing  values 
of  d/r.   The  first  step  in  this  investigation  could  be 
further  experimentation  with  the  configuration  described  in 
this  report  at  various  Reynolds  numbers.   Additional  work 
should  also  investigate  the  effect  of  the  plate  length  on 
flow  patterns.   As  the  wind  tunnel  used  in  this  work  has  a 
capability  for  flow  visualization,  a  series  of  experiments  to 
determine  flow  patterns  visually  would  prove  to  be  valuable. 
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APPENDIX  A 
HEAT  TRANSFER  AND  PRESSURE  DATA  REDUCTION 

The  equations  and  procedures  used  to  reduce  the  heat 
transfer  and  pressure  data  are  shown  below.   The  results  of 
each  experiment  are  presented  on  polar  plots  of  Nu  and  Cp 
versus  Theta  in  section  V. 

A.    HEAT  TRANSFER  EXPERIMENT 

The  relationship  between  power  dissipated  (by  convection 
and  radiation  heat  transfer),  and  power  input  by  Joulean 
heating  is  the  same  as  developed  by  Field  (2)  and  used  by 
Field  and  McComas  (1). 


t  '  RA(T-T  )(0.293) 

RT 1 1 

where:    h.  =  total  heat  transfer  coefficient, ~ 

*  hr  ft   °F 

V  =  voltage  appl ied, vol ts 

R  =  total  resistance  of  Temsheet,  Ohms 

2 
A  =  total  heated  surface  area  of  Temsheet,  Ft 

T   =  surface  temperature  of  the  cylinder  as 
determined  by  liquid  crystals,  °F 

T   =  air  stream  temperature,  °F 

no  * 


0  293  wa-tt-hr  -   746  watts 


BTU     2545  BTU/hr 
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The  Nusselt  number  is  defined  as: 


Nu  = 


hcd 


where 


Nu 

hc 

d 

k 


=  Nusselt  number 

=  convection  heat  transfer  coefficient, 

=  cylinder  diameter, ft 

=  0.015  hrB^oF   (air  at  75°) 


BTU 


hr  ft2oF 


Since  Nu  is  defined  by  the  convective  heat  transfer 
coefficient,  it  is  necessary  to  subtract  the  radiation  effects 
from  the  total  heat  transfer  coefficient 


hc  ■  ht  -  hr 


This  correction  is  made  using  the  equivalent  radiation  heat 

transfer  coefficient  h.  where  h   is  defined  as: 

r         r 

h   -?;.2o(T  +  Tj  (T2  *    T2) 


where 


equivalent  radiation  heat  transfer 

...  •   +    BTU 
coefficient  — 


=  0.1714  x  10 


hr  ft2°F 
-8    BTU 


hr  ft  °R 


2 — -t   (Stef  an-Bol  tzman  Const) 


1-2 


0.9  radiation  exchange  factor  between  the 
cylinder  and  the  surroundings 

cylinder  surface  temperature.  °R 

air  stream  temperature, °R 
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A  sample  calculation  for  data  obtained  at  a  Re  =  153,000 
and  d/r  =  0.50  is  shown  below: 


6=68  degrees,  V  =  24.65  volts, T  =  111.0°F   T   =  63.5°F 

h  =  


24. 652   volt2 


(13.2rt)(0.314:ft2)(lll-63.5oF)(.293wat^T,h,r) 


BTU 


ht  =  10.5 


BTU 
hr  ft2oF 


h   =  (0. 1714x1  0"8  ^~ — 3-)  (  .  9)  (  571  +  523. 5°R) 

r  hr  ft^°RH 

(571 2  +  523. 52R2) 


hr  =  1.0 


BTU 


hr  ftT°F 


ht"hr 


9.5 


BTU 


hr  ft2oF 


BTU 


Nu  = 


(9.5    DIU9   )(.25  ft) 
hr  f-T°F 

(.0150  ■  BT"C) 


hr  ft°F 


u  =  158 


B.    PRESSURE  COEFFICIENT  EXPERIMENT 

As  described  previously,  the  pressure  coefficient  is 
defined  as: 

P-P      P-p 


Pi,   ii2   P  -p 

h   p  U     s   - 
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where : 


C   =  pressure  coefficient 

P 

3 

p   =  density  of  air,  lb/ft 

U^  =  air  stream  vel oci ty, ft/sec 

P   =  pressure  on  cylinder  surface,  in 


H20 


P   =  stagnation  pressure,  in 
P   =  air  stream  pressure,  in 


H„0 


Ho0 


A  sample  calculation  for  data  obtained  at  a  Re  =  153,000 
and  d/r  =  0.50  is  shown  below: 


e  =  70°   P  =  10.0   P   =  6.0   P   =  1 .86 

°°         s 


10.1  -  6.0 
1  .86  -  6.0 


C   =  -  0.99 
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APPENDIX  B 
UNCERTAINTY  ANALYSIS 


Because  of  the  uncertainties  in  the  measured  variables, 
it  was  necessary  to  perform  an  uncertainty  analysis  in 
order  to  ascertain  the  degree  of  uncertainty  in  the  experi- 
mental results.   The  method  of  Kline  and  McClintock  (7) 
was  used. 

The  following  estimated  uncertainties  existed: 
vol tage ,  V  +  0. 1  vol  ts 
resistance,  R  +  0.2  ohms   (20 
temperature  AT  +  1.0°F 


area 


A  +  0.01  ft' 


(20 

1  ) 

(20 

1  ) 

(20 

1  ) 

(20 

1  ) 

(20 

1) 

(20 

1  ) 

(20 

1  ) 

(20 

1  ) 

diameter    d  +  0.001  ft   (20 

radiation       ^ 
exchange  factor  (r  +  0.05 

wall  manometer   H  +  0.05  in 

manometer  bank    S  +  0.05  in 

Total  Heat  Transfer  Coefficient,  h. : 

As  shown  in  Appendix  A,  the  equation  used  to  calculate 
the  total  heat  transfer  coefficient  was: 

V2 


"t  -  RA(T-T  ) 

(since  the  conversion  from  watts  to  BTUs  is  a  constant 
it  does  not  enter  into  the  uncertainty  analysis). 
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For  data  obtained  at  a  Reynolds  number  of  153,000 
and  d/r  =  0.50: 


24. 6^  +  0.1  volts 
13.2  +  0.3  ohms 
47.5  +  1.0  °F 


V 
R 
AT 
A  =  0.315  +  0.001  ft 


AT 


h. 


-r--=0.03   or  u,   =  +  0.3 
ht  ht   " 


ht  =  10.5  +  0.3 


In  some  regions  v/here  T  =  +_  2.0  °F,  then: 

~fc=  0.05   or  u.  =  0.5 
h*  h 

h  =  10.5  +  0.5 

All  measured  quantities  were  read  with  the  uncertainties 
shown  above,  with  the  result  that  the  calculated  value  of 
total  heat  transfer  coefficient  had  a  relatively  constant 
value  around  the  circumference  of  the  cylinder. 
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Since  the  angular  location  of  the  color  bands  does  not 
enter  into  the  calculation  of  H.,  the  angular  uncertainty 
has  not  been  included  here.   This  uncertainty  is  estimated 
to  be  +  3  degrees . 
Equivalent  Radiation  Heat  Transfer  Coefficient,  h  : 

The  equation  used  to  calculate  H   was: 
^  r 

h   =  o(P(1   +  T  )(T2  +  T2) 

Y*  oo  oo 


defining      T      =    1/2    (T   +   T    ) 
3        m  °° 


hr   =    4o*  T; 


since  o    is  a  constant,  it  does  not  enter  into  the  uncertain 
ty  analysis. 

For  the  same  data  used  to  perform  the  uncertainty 

analysis  for  h,,  the  uncertainty  in  h   was  found  to  be: 


m 


(0 


-r-1-  =  0.06  or  w.   =  0.06 
hr  hr 


h   =  1.00^  0.06 
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Nusselt  Number,  Nu : 

The  equation  used  to  calculate  Nu  was 


Nu  = 


h   d 

c 


where  h  =  h.  -  h 
c    t    r 


since  k  is  a  property  of  air,  and  its  uncertainty  is 
negl i  gi  bl e  : 


'Nu 


Nu 


CO. 


2      2 

CO.     +   CO, 

ht    hr 


4   > 


CO 


Nu 


Nu  =  °'03   or  wNu  =  t   5 


Nu  =158+5 


While  this  analysis  has  shown  that 


'Nu 
Nu 


=  0.03,  and 


e  =  +  3  degrees,  it  should  be  noted  that  this  may  be 
misleading.   Because  there  are  portions  of  the  cylinder  which 
exhibit  different  flow  phenomena,  there  are  different  types 
of  color  bands  exhibited.   The  isothermal  region  has  no 
distinct  bands,  and  the  mottled  colors  which  appear  cover 
a  wide  area.   The  liquid  crystals  were  calibrated  with 
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respect  to  the  onset  of  each  color.   The  result  is  that 
for  areas  of  pulsating  or  mottled  colors,  the  uncertainty 
of  surface  temperature  and  location  are  much  higher  then 
forregions  in  which  a  distinct  color  gradient  exists. 


Reynolds  Number,  Re : 

The  equation  used  to  calculate  Re  was: 

U   d 
Re  =  


The  value  of  U   is  directly  proportional  to  the  sqaure 
root  of  the  wall  manometer  reading,  therefore: 


10 


U 


u 


=  0.5 


'H 


H 


since  v  is  a  property  of  air,  and  its  uncertainty  is 
negligible: 


Re 


0.5  uu\  2 

nr-1 


u 


Re 


=  0.011   or  wD   =  1700  «  2000 
Re  Re 


Re  =  153,000  +  2000 
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Pressure  Coefficient,  C  : 

The  equation  used  to  calculate  C   was 

P 


P-P 


P  -P 

«  s 


A  P 
A~P~ 


both  pressure  differences  were  measured  on  the  same  manometer 
bank.   For  Re  =  153,000  and  d/r  =  0.50: 


P    =    10.1         P      =    6.0        P      =    1 .86 

°°  s 


■ap 


■ap. 


ap 


AP, 


(O, 


0.017      or   u>c      =    0.017    «  0.02 
P 


C      =-0.99    +    0.02 
P 
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APPENDIX  C 
HOT  WIRE  ANEMOMETER  CALIBRATION 

Prior  to  any  experimental  runs,  the  TSI  hot  wire 
anemometer  was  calibrated  using  a  controlled  air  stream  and 
a  precision  manometer.   In  order  to  assure  accurate  calibra- 
tion, the  probe  was  mounted  in  the  wind  tunnel  traversing 
mechanism  and  the  calibration  was  done  in  the  wind  tunnel 
test  section.   In  this  way  all  leads  and  connections  used 
for  experimental  runs  would  be  in  the  circuit  for  calibration 

The  calibration  apparatus  provided  a  controlled  air 
flow  whose  velocity  was  calculated  from  the  pressure 
differential  provided  by  an  al cohol -f i 1 1 ed  manometer.   The 
velocity  calculated  was  plotted  against  the  voltage  indicated 
by  the  hot  wire  anomometer  unit  (Figure  20).   Calibration 
was  done  by  increasing  the  calibration  apparatus  air  stream 
velocity  so  that  the  alcohol  in  the  manometer  would  indicate 
even  values  (in  mm).   The  manometer  inclination  and  readings 
were  recorded,  along  with  the  voltage  indicated  by  the  hot 
wire  anemometer. 

The  equations  used,  data  taken  and  calibration  curves 
are  shown  below.   Calibration  checks  were  made  periodically 
throughout  the  course  of  this  work,  and  corrections  made  to 
account  for  the  aging  of  the  probe. 
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From  Bernoulli's  Equation 


V 


2  g  AP 


where : 


U 

00 

P 

AP 


air  stream  velocity  ft/sec 

32.17  ft/sec2 

3 

air  density  lb/ft 

pressure  differential  lb/ft' 


The  pressure  differential  from  the  manometer  was  calculated 
from  the  following  relations: 


%0  g 
P  =  — = h 


and 
where 


0.789 
0.80 


ch' 


0.789 
0.80 

c 
h' 


ratio  of  specific  gravity  of  manometer 
fluid  used,  to  specific  gravity  of 
fluid  intended  for  manometer 

manometer  inclination  factor 

height  of  miniscus  in  inclined  manometer  mm 


hence 


where : 


U      =    13.24    /h      ft/sec 


h      =   mm   of   al cohol 
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Turbulence  intensity  was  determined  by: 


where 


V 
e 
V, 


4  V  e 

2     2 
V2  -  V2 


turbulence  intensity 

voltage  indicated  by  hot  wire  anemometer 

rms  voltage  variation 

extrapolated  value  of  voltage  corresponding 
to  zero  velocity 
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Voltage  vs  Velocity 


Figure  20.   Voltage  sqaured  versus  the  square  root  of  velocity 
for  hot  wire  anemometer  calibration. 
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APPENDIX  D 
WIND  TUNNEL  VELOCITY  PROFILE  AND  TURBULENCE  INTENSITY 

In  order  to  assure  that  test  results  would  be  accurate, 
it  was  necessary  to  show  that  the  velocity  profile  in  the 
test  section  was  constant  and  that  turbulence  was  small 
enough  to  have  no  effect  on  test  results. 

The  hot  wire  anemometer  was  used  to  probe  the  test 
section  at  various  locations.   The  test  section  coordinate 
axis  is  shown  in  Figure  21.   The  test  cylinder  was  located 
at  X  =  13.71,  Y  =  -2.33.   Velocity  profiles  were  taken  at 
the  following  locations: 

Test  section  entrance:   X  =  17.75  -9 .O^Y^+9 .0   Z  =  0.0 

X  =  17.75  Y  =  -2.2£+4.0£Z  =4.0 
X  =  -1 .71  -3.02  3<_Y<  +  9.0  Z  =  0.0 
X  =  -13.71  -0.825<Y<+9.0  Z  =  0.0 


Plate  leading  edge: 


Top  of  the  cylinder: 
(plate  attached) 

Top  and  bottom  of  cyl .  :  X  =  -13.71  -10.0<Y<  +  9.0  Z  =  0.0 
(plate  at  d/r  =  5.33) 

All  velocity  profiles  were  taken  at  a  tunnel  speed  of 
98  ft/sec. 

As  shown  in  Figures  22  and  23  the  air  speed  increased 
from  10  to  30  percent  on  either  side  of  the  cylinder  in  the 
case  of  the  plate  at  d/r  =  5.33,  and  a  like  amount  over  the 
top  of  the  cylinder  in  the  case  of  the  attached  plate.   The 
velocity  at  the  test  section  entrance  was  essentially  constant 


64 


Turbulence  intensity  was  determined  at  the  test  section 
entrance : 

X  =  17.7  -S.0<y<+5.0    1    =    0.0 

for  three  Reynolds  numbers;  75,000,  100,000,  and  153,000. 
The  average  value  of  turbulence  intensity  for  Re  =  153,000 
was  0.00239,  or  about  0.2%.  The  turbulence  for  the  other 
Reynolds  numbers  was  slightly  higher,  but  of  the  same  order 
of  magnitude  (see  data).   The  value  of  turbulence  intensity 
was  calculated  in  the  same  manner  as  for  the  hot  wire 
anemometer  probe  calibration  described  in  Appendix  C. 
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FLOW 

DIRECTION 


TEST    SECTION    TRANSLATOR 
POORDINATE     AXIS 
Figure    21.       Isometric    of    the    test    section    coordinate    axis 
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